Flexible Linear Polyelectrolytes in Multivalent Salt Solutions: Solubility Conditions 
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Polyelectrolytes such as single and double stranded DNA and many synthetic polymers undergo 
two structural transitions upon increasing the concentration of multivalent salt or molecules. First, 
the expanded-stretched chains in low monovalent salt solutions collapse into nearly neutral compact 
structures when the density of multivalent salt approaches that of the monomers. With further ad- 
dition of multivalent salt the chains redissolve acquiring expanded-coiled conformations. We study 
the redissolution transition using a two state model [F. Solis and M. Olvera de la Cruz, J. Chem. 
Phys. 112 (2000) 2030]. The redissolution occurs when there is a high degree of screening of 
the electrostatic interactions between monomers, thus reducing the energy of the expanded state. 
The transition is determined by the chemical potential of the multivalent ions in the solution /j, 
and the inverse screening length k. The transition point also depends on the charge distribution 
along the chain but is almost independent of the molecular weight and degree of flexibility of the 
polyelectrolytes. We generate a diagram of fj, versus k 2 where we find two regions of expanded con- 
formations, one with charged chains and other with overcharged (inverted charge) chains, separated 
by a collapsed nearly neutral conformation region. The collapse and redissolution transitions occur 
when the trajectory of the properties of the salt crosses the boundaries between these regions. We 
find that in most cases the redissolution occurs within the same expanded branch from which the 
chain precipitates. 

I. INTRODUCTION 



The precipitation and dissolution of linear polyelec- 
trolytes with the addition of multivalent salt or particles 
have been extensively studiedtflj. In DNA the precipi- 
tation provides a promising mechanism to "pack" long 
DNA, a critical problem in gene therapy. Moreover, it 
is correlated with highly accelerated rates of DNA re- 
naturation and cyclisatioril3'tll. Therefore, it is impor- 
tant to determine the concentration of multivalent salt 
or particles at which the precipitated DNA dissolves in 
the solution. Furthermore, it is crucial to determine the 
effective charge of the DNA in the precipitated and dis- 
solved states to determine their interaction with cells. 

It has been shown experimentally that the transi- 
tions are nearly independent of the details of the lin- 
ear polyelectrolyte (such as charge density, degree 
flexibility, molecular weight, structure function, etc.) 
That is, flexible (single stranded DNA and Polystyrene 
Sulphonate) and semiflexible chains (double stranded 
DNA) can be described by the same thermodynamic 
model even though flexible chains collapse into amor- 
phous dense spheres—and long semiflexible chains into 
toroid conformationsE3. 

We recently develop a thermodynamic model that ex- 
plains the universal nature of the multivalent counterion 
induced precipitation transition in low concentration of 
monovalent saltaJ. In this paper we analyze the redis- 
solution transition by extending our two state model to 
include large concentrations of multivalent salt. 
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FIG. 1. Location of the collapse and redissolution transi- 
tions in a diagram of the logarithm of the concentrations of 
monomers <j> an d of multivalent salt m, for dilute, strongly 
charged polyelectrolytes. Lines LI and L2 form the bound- 
ary between the expanded-stretched state and the collapsed 
state, and line L3 is the boundary between collapsed and 
expanded-coiled states. The line P is a guide used in the 
text and represents the increase of multivalent salt at fixed 
monomer concentration. 



The experimental diagrams of the logarithm of con- 
centration of monomers <p versus the logarithm of the 
concentration of multivalent salt m have three regions 
(A, B ajid C) separated by 2 transitions lines (L1-L2 
and L3)BB, shown schematically in Fig. 1. Region A 
corresponds to polyelectrolytes dissolved in water with 
expanded-stretched conformation at low m values. Re- 
gion B at intermediate multivalent salt concentrations 
m is a solution of collapsed polyelectrolytes (sometimes 
considered as a coexistence of two phases: one rich and 
one poor in polyelectrolytes). Region C at large m val- 
ues contains polymers dissolved in water in expanded- 
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coiled conformations. The transition line between the 
A and B regions, L2, at low monovalent salt concentra- 
tions is linear with a slope comparable to X/z, where z 
is the valence of the counterions of the added salt. In 
the presence of monovalent salt there is a nearly hori- 
zontal transition line, LI, at very low concentrations of 
monomers. The transition line between regions B and 
C is horizontal (i.e., independent of the concentration oi 
monomers) in the full regime. In our previous paperQ 
we discussed the transition between regions A and B and 
explained the transition in the regime L2 as the creation 
of collapsed conformations of polymers whose charge is 
compensated by condensed multivalent salt. Our model 
also predicts the transition LI in the presence of large 
amounts of monovalent salt. In this paper we study the 
transition between the regions B and C in solutions with 
very low (negligible) monovalent salt concentrations. 

It is well documented that in polyelectrolyte solutions 
a fraction of the counterions condensed along the poly- 
electrolytes to decrease their electrostatic energy^. This 
ion condensation is crucial to understand the precipi- 
tation and redissolution transitions. In this paper we 
compute the fraction of condensed ions as the multi- 
valent salt concentration increases in the two possible 
polyelectrolyte conformations: collapsed and expanded. 
Our two state model is based on the fact that in salt 
free and/or low ionic strength solutions, ion condensa- 
tion leads to two possible conformation in linear flex- 
ible polyelelctrolytescJ: expanded-stretched with a re- 
duced effective charge and collapsed jaearly neutral with 
the cohesive energy of an ionic glasso. Indeed, multiva- 
lent counterions lead to large cohesive energies for nearly 
neutral collapsed structures in region B explaining the 
transition L2. Since the magnitude of the electrostatic 
energy of a collapsed chain in region B is much larger 
than any entropic energy increase (due to the decrease 
of degrees of freedom resulting from the compaction), 
our model also describes the precipitation of semifiexible 
and rigid-rod polyelectrolytes into toroids and bundles, 
respectively, explaining the universality of thej-precipita- 
tion transition discussed in our previous workQ. 

At the precipitation transition we predict that the ex- 
panded chains are slightly charged, and with further ad- 
dition of multivalent salt, the collapsed chains (in region 
B) become practically neutral, in excellent agreement 
with recent electrophoresis experiments!! As a result, the 
aggregates interactions in the dilute regime are negligi- 
ble justifying a monomolecular collapse model to describe 
also multi-molecular aggregation. 

The effective charge and size of the expanded con- 
formation (in region C) changes with further addition 
of multivalent salt. The transition L3 occurs at high 
salt concentrations, and in such conditions the chains 
are expected to obey random walk or self-avoiding 
walks statisiics due to the screening of the electrostatic 
interactionsti Screening is responsibkj-pf the redissolu- 
tion, as suggested in various modelsQ'Elcl. We show here 
that the redissolution is actually very sensitive to the re- 



lation between the chemical potential of the condensing 
multivalent particles in the solution \i and the inverse 
screening length n, while nearly insensitive to the de- 
gree of flexibility and molecular weight of the polyelec- 
trolytes. Screening, for example, may be reduced due to 
multivalent-monovalent ion associating in the solution, 
strongly affecting the effective charge of the expanded 
chains. We determine the effective charge of the chains 
as a function of fi and k including the finite size of the 
ions and the discrete nature of the charge distribution 
along the polyelectrolyte, which is essential to compute 
the electrostatic energies of both the collapsed and ex- 
panded statesQ. In a /i - k diagram the collapsed region 
lies between two branches of expanded states: one ex- 
panded with a reduced effective charge and the other 
expanded with an inverted charge. In most common 
situations the chains precipitate and redissolve within 
the same branch; i.e., with a reduced effective charge 
of the same sign of the bare charge. Our model predicts 
a redissolution transition independent of polyelectrolyte 
concentration and describes the re-dissolution of flexible, 
semifiexible and rigid polyelectrolytes. Indeed, the col- 
lapsed and expanded states discussed here are akin to the 
multi-molecular precipitated and dissolved states, respec- 
tively, observed in many polyelectrolytes^, explaining the 
rather universal form of the diagram, Fig. 1. 

In the Section 2 we summarize the assumptions of our 
work, and provide a description of each of the phases 
of the system. In section 3 we construct the free ener- 
gies of two types of representative states and in Section 
4 we give the result of the free energy minimization. In 
Section 5 we summarize our results, compare them with 
experiments and emphasize the elements of the theory 
that have not yet been subject to experimental results. 
We end with a brief conclusion in section 6. 



II. THEORETICAL MODEL 

As in our previous workEl, we rely in a two state de- 
scription of the system. Instead of considering all pos- 
sible states of the system (describing them, for example 
by a changing scaling exponent), we simply assume that 
the minimum lies in one of two extremes: collapsed or 
expanded conformations. The expanded conformation, 
however, is a function of the concentration of salt in the 
solution. For example, at low monovalent salt concentra- 
tions the chains are expanded rods, and at large concen- 
trations they are expanded coilstl. The free counterions 
and solvent can be integrated out provided we take into 
account the electrostatic interactions between the non- 
condensed counterions in the solution. 

The condensed (collapsed) state is treated as an ionic 
glassQ so that its energy can be approximately calculated 
by techniques from solid state physicsEH; i.e, we assume 
that fluctuations play no important role and can be ne- 
glected. For a very large chain, at low enough tempera- 



2 



tures, the bulk of the collapsed state should acquire an 
almost crystalline structure. The finite size of the chain, 
the connectivity of the monomers, and finite tempera- 
ture effects must induce defects into the structure and 
the local structure should be glass-like, where the perfect 
order has not been achieved. An ionic glass structure was 
indeed observed in recent numerical studies of polyclcc- 
trolytes in bad solventsEEl. 

In both states of the systems, the collapsed and ex- 
panded states, it is necessary to keep track of the finite 
size of the counterions and monomers. This is clear in 
the case of the collapsed state, but even in the expanded 
states there is an important contribution arising from in- 
teractions with the nearest neighbors, and this energy 
can only be.-calculated by considering the finite size of 
the particles!! For example, a naive coarse-grained model 
that takes into account only effective charges would as- 
sign a zero energy to a neutral cluster, while it can ac- 
tually have a very strong cohesive energy. The redissolu- 
tion into an expanded state occurs at high salt concen- 
trations and the effective interactions between monomers 
are screened. In these conditions the energy associated 
with interactions between far away segments of the chain 
can be treated as a small perturbation with respect to 
that of the neutral polymer state. 

The properties of the salt at different concentrations 
are summarized into two parameters: an effective chem- 
ical potential fi, and an inverse screening length n. We 
will present our results for general values of these param- 
eters and in particular consider the case were the param- 
eters are related by a Debye-Huckel law. We also discuss 
the effect of multivalent-monovalent ion association, the 
Bjerrum modelES, in the diagram. 

Let us consider the series of states of a single polymer 
chain along a line of constant polymer concentration (line 
P in Fig. 1) with increasing multivalent salt concentra- 
tion. We start in a state in region A with few counterions 
condensed. The condensed counterions lie very close to 
the chain, and cannot be considered as mixed in the sol- 
vent. The effective charge induces repulsions between 
the segments of the chain, and the chain acquires extra 
stiffness, making it into a rod-like structure. 

When a small quantity of multivalent salt is added to 
the system, it dissociates, and the multivalent counteri- 
ons are absorbed to the chain because the electrostatic 
energy to be gained from the condensation is larger (per 
unit charge) than that of the monovalent counterions. 
Furthermore, the entropic loss due to condensation is 
lower for the multivalent counterions. When the charge 
of the added salt roughly equals the bare charge of the 
polymers the effective charge of the chains quickly goes 
to zero. At this point most charges of the chain are com- 
pensated by the condensed multivalent counterions and 
a few monovalent ones. Since a group of charges with 
total charge zero has a minimal energy when they are 
set in a compact configuration, a collapsed conformation 
is acquired. This transition is of course also possible for 
the case of pure monovalent salt, but it requires cither 



larger concentrations or lower temperatures to offset the 
entropic contributions. In other words, the determining 
factor for the transition is the strength of the electrostatic 
interaction, and the introduction of multivalent salt acts 
as a jump in parameter space and is not equivalent to a 
smooth increase in monovalent saltH. 

In the collapsed state the absence of uncompensated 
charges in the polymer reduces its solubility in water 
and polar solvents. It is not clear from experiments, 
however, whether the collapsed flexible polymers aggre- 
gate into large structures or if they form a dense solu- 
tion of individually collapsed chains. Though we assume 
monomolecular collapse, the predicted transition points 
should not exhibit strong dependence on the number 
of chains aggregated since if the collapsed state is neu- 
tral, the residual interactions between chains (promoting 
or opposing multimolecular aggregates) are energetically 
small. 

Since the collapsed structure is almost neutral, the in- 
crease in salt concentrations has little effect in the ab- 
solute value of the energy of the chain. The change in 
the environment affects only the interaction between the 
positive or negative small excess (or defect) of charge in 
the polymer, which would be located at the surface of the 
collapsed structure. The energy of the charges packed in- 
side the structure is dominated by the interaction with 
their nearest neighbors and is not subject to screening 
(there are no floating ions between them). On the other 
hand, the energy of an expanded chain with non-zero ef- 
fective charge will be greatly reduced by the increased 
screening. 

For the purposes of comparison with the collapsed 
structures, we can consider a very rigid chain that cannot 
be collapsed into a compact structure. We can ignore for 
the moment the fact that this chain would likely form 
a bundle with other chains, and we can assume that we 
have isolated it by perhaps mechanical means. At higher 
multivalent salt concentrations, the rigid polymer contin- 
ues absorbing ions, and thus it can become neutral and 
even overcharged (effective charge with opposite .sign-to 
the bare charge, as it occur in many systemso'OcJo) . 
With further increase in the amount of salt the trend can 
be reversed, and the chain can again become, first neu- 
tral, and then simply charged (same charge sign as the 
bare charge). The increase of salt concentration reduces 
the entropic penalty for the condensation of counterions. 
The counterions continue condensing into the polymer 
as long as there is a reduction in electrostatic energy due 
to the condensation. The increase in condensation stops 
when the typical screening length of the free counteri- 
ons is of the order of the monomer separation, the free 
counterions can reduce their energy by associating with 
each other, and a fraction of them will do so instead of 
associating with the polymer. 

At large salt concentrations, the chain reverts from col- 
lapsed to expanded, giving the transition B — C. In the 
expanded state the chain needs not to be neutral. The 
screening of electrostatic interaction under these condi- 
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tions strongly reduces the interaction between the un- 
compensated charge of the polymer, and it is clear that 
the repulsions will not be strong enough to create a rod- 
like state. In these conditions an expanded coil state is 
more likely. In terms of the effective charge of the poly- 
mer we find that two scenarios are in principle possible. 
If the chemical potential of the counterions in the sol- 
vent is low enough while the screening length is reduced, 
the expanded state is overcharged since there is an effec- 
tive decrease of energy by increasing the condensation to 
the chain, as recently discussed by Nguyen et al.Q. On 
the other hand, when the screening length is small, it is 
possible that the interactions between coions and coun- 
terions in the solvent become important, and it is much 
harder for the chain to acquire the counterions. This is, 
the chemical potential for the extraction of counterions 
from the solvent is large and therefore the effective charge 
at the redissolution point has the same sign as the bare 
charge. We show in the next section that though there 
is a possibility for overcharged polyelectrolytes at the re- 
dissolution, in the sense that this state can be preferred 
to a collapsed one, the salt properties required for this 
situation (/i, k) are unlikely to be realized in practice for 
water soluble multivalent ions. 



III. CALCULATION OF FREE ENERGIES 

We assumed that the polymers are monodispcrse with 
monomer number N. Each monomer is charged posi- 
tively, and is monovalent. The original counterions of 
the polymer also have valence 1. The polymer density is 
<f>. We measure all distances with respect to the monomer 
size b (so that b — 1 in our units). It is necessary for more 
precise calculations to also give more detailed geometri- 
cal information such as the aspect ratio of the monomers 
and the respective size of all the ions. Here we only need 
to further consider, besides the size of the monomer, the 
radius of the multivalent counterions a. The monovalent 
salt concentration s is set equal to zero, but there are 
monovalent counterions originally dissociated from the 
polymer with density (j), and the multivalent salt, with 
valences z : 1 has concentration to. All energies are mea- 
sured in units of fcgT, (the Boltzmann constant times the 
temperature), and consider only energies per monomer. 
We form the dimensionlcss Bjerrum number B, equiva- 
lent to the Manning parameter, as 



4£ vr6fc B T' v ' 

where £o is the permitivity of the solvent (water) . The 
typical value of this number at room temperature in 
strongly charged polyelectrolytes is larger than 1 and as 
high as 4.2 for double stranded DNA. We assume here for 
simplicity that B is a constant, independent of the num- 
ber of locally associated monomers and ions. Given that 



in both states water is always present, the corrections in 
the values of e are expected to be minora. 

In both expanded and collapsed states we describe the 
system by a particular geometry and a given fraction x 
of condensed charge. The free energy is separated into a 
contribution from the solvent and free ions, and another 
part from the interaction of the monomers and charges 
condensed to them. The fraction condensed is obtained 
self-consistently from minimization over all possible con- 
densed fractions. The form of the entropic term is the 
same in both cases, but it is easier to consider the elec- 
trostatic contributions separately. 

A. Free counterions contributions 

It is convenient to separate the region occupied by the 
polymer from the rest of the solvent. To calculate the 
free energy of the solvent region, we first calculate the 
entropic contribution. We assume that the charge con- 
densed to the chain — xN is composed by x s N monova- 
lent ions and x m N/z multivalent ions. Thus there are 
{s/<j>) + \ — x s free monovalent counterions per monomer, 
and thus the energy per monomer in fc^T units is 

F.i = - — ) ln(m - x m ^) + (i + 1 - xs) ln(s + (1 - x s )<t>). 
q> z z <p 

(2) 

In this expression we omitted the contribution of the ions 
with the same charge of the monomers, since it is as- 
sumed that they do not condense, and thus have a con- 
stant contribution. We have also omitted the contribu- 
tion of the solvent, since the concentration of the salts is 
small compared to the solvent. It is possible to add in 
this expression an interaction term between the solvent 
and the ions, but for simplicity we neglect these correc- 
tions here; its effect would be only to shift the position 
of the L3 line in the phase diagram. It is important to 
consider, however, the effective free energy of interaction 
between all free ions. A simple way to do this is to use 
the modified Debye-Huckel form: 

F S 2 = -. — - ( ln(l + na) - na + -(ka) 2 ) . (3) 

47T0 \ 2 J 

In this expression a is the hard-core radius of the salt 
ions and k is the inverse screening length produced by 
the salt: 

K 2 = 4irB(z 2 (m - x m (j)/z) + zm + (2s + (1 - x s )4>)). 

(4) 

In the very dilute polymer limit, we can simplify these 
expressions, and use instead 

k 2 = Airz(z + l)Bm, (5) 

and expand the total free energy of the solvent F s = 
F s i + F S 2 in a scries in the overcharge fraction y = x — 1. 
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The constant term is irrelevant and we retain only the 
linear term: 

F s = (d y F s )y = fiy, (6) 

where p is now the effective chemical potential for the 
extraction of ions from the solvent into the condensed 
region. Using the same high salt limit as before we ob- 
tain: 

1 1 K 3 

(i = In m H r . 

z 8t:{z + l)m 1 + Ka 



B. Expanded state energy 



In our previous worku, we considered only the energy 
of rod-like conformations. Since our interest is the regime 
with large amounts of salt, the screening strongly re- 
duces the rigidity of the chain due to excess charge. It is 
then adequate to consider a Gaussian conformation (self- 
avoiding walk conformations give very similar results, as 
it can be seen from the functional form of the free energy 
computed below). The electrostatic contribution to the 
free energy is calculated as a function of the condensed 
charge, and we set the neutral state x = 1 as reference. 
The free energy is expressed as the sum of three terms. 
A constant term reflects the energy of the neutral state, 
a term linear in the excess charge y — (x — 1) appears 
due to addition or extraction of a charge into or from 
the neutral state. Finally, the interaction between seg- 
ments of the chain with excess charge gives rise to a term 
quadratic in y. We write this as: 

F e e ^9 e Q +g e 1 y + 9 c 2 y 2 . (7) 

A coarse-grained model clearly gives rise to the last term, 
and the first two can be understood, in that context, as 
regularizations that take into account the finite size of the 
particles involved. In the limit of large amounts of salt- 
added, most of the condensed charge comes from multi- 
valent ions, and thus, we will take x = x m , and x s = 0. 
The terms are graphically represented in Fig. 2. 

For the calculation of the energy of .the neutral state we 
basically apply the Wigner principleO: in a dense ionic 
system, the free energy can be approximated by the con- 
tribution of nearest neighbors that effectively cancel the 
charge of the particle. In the case of the collapsed state, 
below, it is more suitable to consider an approximation 
based on the calculation related to a full infinite lattice. 




FIG. 2. Graphical representation of the terms required to 
calculate the energy of the system in both states. The go term 
measures the energy per monomer of the chain in a neutral 
state. i?i measures the energy of adding an additional charge 
to the neutral conformation. <?2 is the energy of interaction 
between added charges and is obtained by a coarse graining 
procedure. All the excess charge in the collapsed state (b) is 
assumed to lie on its surface. 

The basic neutral cluster of the linear polyelectrolyte 
consists of one multivalent counterion, and z monomers. 
Since all particles in the cluster (with charges qi, both 
monomers at distances r.y from each other, all roughly 
equal to the monomer size r = b = 1, the electrostatic 
energy of the cluster is 

E = Y / B^ = -B(z 2 + z(z-l)/2) (8) 
&i 3 

Dividing by z monomers, we obtain: 

g e = -B(z+l)/2. (9) 

When an extra multivalent ion is added to this cluster, 
the counterion can have roughly the same interaction 
with the chain monomers, but can locate itself away of 
the first counterion, say, by putting a monomer in be- 
tween them so that r rj 2. Thus we obtain a simple 
estimate for the energy of the new charged cluster as: 

E = B(-2z 2 + z(z -l)/2 + z 2 /2), (10) 

and 

9\ = -Bz/2. (11) 

Finally, the two-body term can be calculated assum- 
ing a uniform distribution of charge along the Gaussian 
chain. A Gaussian chain has a monomer distribution 
given, approximately, by p(r) = l/(2-7rr), inside the vol- 
ume V limited by the radius of gyration R; this is, for 
r < R = N 1 / 2 . A charge at the center of this distribu- 
tions interacts with other charges via the screened poten- 
tial Be~ Kr jr. The energy for one such charge is given by 
the average over the volume V of this potential: 

r p—Kr f> 

E=(V)= drp(r)B = -(1 - e""*), (12) 

Jv r k 
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and by our previous definition 
„ 1 



.92 



(V). 



(13) 



Notice that while this calculation depends on the as- 
sumed Gaussian conformation, a simple estimate can also 
be done by using a general expansion coefficient for the 
radius of gyration R w N v , (v = 1 for rod-like, v = 2/3 
for self-avoiding walks, and v = 1/2 for random walks). 
In such case, we obtain a scaling for the quadratic preco- 
efficient of the form g\ k^ 1 /"). While these different 
models give rise to different estimates, the transition we 
are interested in occurs when the screening is of the order 
of the size of the monomers, and so, in our units, K ks 1, 
and the estimates at that point are not so different. 



C. Collapsed state energy 

We repeat a similar approach for the calculation of the 
free energy of the collapsed state which is expressed as: 



c„,2 



K = 9o + 9iV + 92V 



(14) 



Inside the collapsed polymer, we can imagine the envi- 
ronment of a charge resembling that of a ionic crystal. If 
the structure of the crystal is known, the energy per atom 
in that crystal can be given in terms of the Madelung 
constantcJ, 



E 



-Me 2 /r 



-MB. 



(15) 



The Madelung constant for an ionic monovalent salt is 
M = —1.747. There are two important differences in 
our case. We would like to consider multivalent ions, 
which reduce the contribution from same-charge inter- 
actions, and on the other hand the connectivity of the 
monomers forces them to be next-neighbors (not the case 
in ionic crystals), and thus increases the same-charge in- 
teractions. A way to obtain a rough estimate for M is to 
consider a simple geometry in which both monomers and 
counterions lie in columns which in turn form a square 
lattice when a cross section is taken. The result does 
not change dramatically when other possible geometries 
are considered. As a rule of thumb, we can use an ef- 
fective Madelung constant that is larger than unity, but 
not much bigger. In the geometry proposed, the number 
of next neighbors of counterion is Az. Summing the in- 
teraction with Az monomers and 4 equally charged next- 
nearest counterions (r = 26) we obtain the following es- 
timate of the energy per cluster of 1 counterion and z 
monomers: 



E = B{-Az 2 + z(z - 1) + -z(Az -z) + 



lAf 
2~Y 



(16) 



where the interactions with charges outside the cluster 
are weighted by 1/2 to avoid double counting. Dividing 



again by z, we obtain the coefficient of the energy per 
monomer: 



ti = -B{l + -). 



(17) 



In this approximation we have underestimated the 
Madelung constant for z = 1, for which we obtain —1.5 
instead of the exact result —1.747, but it is a reasonable 
approximation and shows that the energy of this neu- 
tral collapsed state is smaller than the neutral expanded 
state. 

In both the collapsed and expanded states, the addi- 
tion and extraction of charges are not completely sym- 
metric processes. It is energetically unfavorable to add 
or extract charges from the collapse bulk state, as shown 
below in the calculation of g%. Both processes are how- 
ever better thought of as occurring at the surface_pf the 
collapsed region. We have also argued previouslyO, that 
while the mobility of charges in this type of conforma- 
tions is not so large, the analogy with a conductor is still 
useful: all charge excess resides in the surface. While the 
local conformation at the surface of the collapsed state 
is not equal to the expanded state, it is simpler to bor- 
row the result for the coefficient related to the addition 
of one charge from that case. In both situations the new 
charge is put in close contact with a neutral cluster, but 
sits slightly closer to the opposite charges. Thus, 



g\ = -Bz/2. 



(18) 



Since we assume that the extra charge is sitting at 
the surface, the interaction energy can be calculated as- 
suming a uniform distribution in a shell surrounding the 
collapsed polymer. We assume this shell to be spherical 
with radius R, and obtain the energy per added charge 
as 



B 



2 R 2{3N{1 + 1/z)/Att) 1 / 3 



(19) 



Here we have considered that the N + N/z charged par- 
ticles are closely packed and that they occupy each a 
volume of b 3 . 



IV. COMPARISON BETWEEN STATES AND 
PHASE DIAGRAM 

After the minimization of the free energies of each state 
with respect to their respective condensed charge vari- 
ables, we can compare each other to find the state of the 
chain for given concentration conditions. This minimiza- 
tion can be done directly from the expressions already 
given but it is easier to present simpler formulas if we 
use a set of approximations described below. Numeri- 
cal results presented later reflect these approximations, 
but we have checked that numerical minimization of the 
full expressions do not lead to important qualitative and 
quantitative changes in the results. 



G 



The free energy of each of the states, for given concen- 
trations of monomers and salt correspond to the minima 
of: 



F e 
F c 



(20) 
(21) 



for the expanded and collapsed states respectively. We 
minimize now with respect to the overcharge fraction 
y. Taking the first derivatives with respect to this quan- 
tity, we obtain the (independent) conditions 



d y F e = ix + g\ + 2g c 2 y e = 0, 
d y F° = n + gl + 2gZy° = Q. 

The solutions for the fractions condensed arc: 



y 
y c 



-(p + g1)/2g e 2 , 
-{ l i + gl)/2g c 2 . 



(22) 
(23) 



(24) 
(25) 



Replacement of this values in the general expressions for 
the free energy gives the minimized free energies for each 
state: 



F e *=F e (y e ) = M + ^ 
F c * = F% c )= M + 3o c - 



(fli+M) 2 
4<?f 

(ff!+M) 2 

4ff 2 c 



(26) 
(27) 



For the moment it is convenient to maintain the ex- 
pressions of the free energies as functions of the chemical 
potential instead of the salt concentration. Of the param- 
eters used to describe the different states, only g 2 has an 
explicit dependence on the concentration (through the 
inverse screening length Also, it is important to note 
that the coefficient g 2 is extremely large, and for the val- 
ues of the chemical potential that we need to explore, it 
turns out that the denominator in the expression for the 
overcharge is much larger than the numerator, and the 
overcharging fraction is very close to 0: 



j/ c «0. 



(28) 



This result is clear since we have found that there is a 
very strong penalty for any overcharging (or undercharg- 
ing) in the collapsed state. The effective charge in a gen- 
eral structure is affected by the environment (namely, 
by fi) but the collapsed structure is only very slightly 
susceptible to it. Furthermore, the interaction between 
the charges in the bulk of the collapsed structure is not 
subject to screening by the free particles outside and its 
contribution remains constant with changing chemical 
potential. Thus, the system splits neatly into collapsed 
neutral structures and a homogeneous mixture of solvent 
and non-condensed counterions. It is then not surprising 
that this state adequately represents what in experiment 
is close to be two separate phases. 

Consider now the minimum for the expanded state. 
The coefficient for charge-charge interactions g 2 varies 



very strongly with the salt concentration. It remains 
large when the screening length is large, since in those 
conditions all uncompensated charges interact strongly. 
On the other hand at large concentrations of salt the 
screening drastically reduces the interaction between seg- 
ments of the chain. 

To obtain an expression for the transition point we sub- 
tract the free energies, and define the energy difference 



AF 



(29) 



The system is in the region B, when AF > 0, in region 
C when AF < 0, and the transition curve is defined by 
AF = 0. Using the fact, discussed above, that the elec- 
trostatic energy of the collapsed state remains essentially 
constant, we can write the equation AF = as 



5o 4g| - +%■ 



(30) 



One useful way to use and interpret this equation is to 
consider it as defining a boundary in the \x — k 2 plane that 
corresponds to the properties of the counterions. Within 
this diagram, a particular form of the relation between 
concentration and screening for a counterion defines a 
path in the diagram parametrized by the salt concentra- 
tion to. Examples of this are presented in Fig. 3. for 
certain polymer properties values, and different types of 
/i — k 2 relations. 

At this point is very easy to explain the observed near 
independence of the transition on the molecular weight 
of the chains. Only the coefficients g 2 carry informa- 
tion on the molecular weight, but as we have seen g 2 
drops out of the final equations, and in g 2 the molecular 
weight provides only an exponentially small correction to 
the dominant term given by a function of the screening 
length. 

As shown in Fig. 3, and as it can be seen from the fact 
that the equations defining the transition are quadratic 
in fi, the expanded state splits into two regions, each 
with non-zero effective charge, but with opposite signs on 
each. This diagram is obtained using the approximations 
presented above, and a more precise determination of the 
lines of transition will involve consideration of the mono- 
valent counterions, the finite amount of monomers, and 
the rod-like state obtained at very low screening. When 
these considerations are taken into account, the trajec- 
tory described by the added salt in the fi — k 2 plane will 
start within the expanded region, cross into the collapsed 
region (giving rise to the L1-L2 line), and continue there 
until crossing again, redissolving, into one of the branches 
with expanded states. The starting point of the trajec- 
tory, with zero multivalent salt added to = 0, is located 
in the branch with natural charge, at a finite k value (due 
to the monovalent counterions), and at \l — oo. Figure 3 
illustrates only the redissolution. 

In Fig. 3 we show three different curves defining the 
properties of the counterions: a pure Debye-Huckel case 
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where the size of the counterions are neglected, a modi- 
fied Debye-Huckel curve with a — 0.2 (that is a = 0.26), 
and a third curve that neglects the interactions between 
the counterions when they are free considers only the en- 
tropic term. The reversal of direction in the fi axis of 
the parametric curve for the Debye-Huckel cases occurs 
because of the onset of important favorable interactions 
between the free charges. Thus it is clear that the tra- 
jectory can end in the naturally charged branch of the 
expanded region since there are more free counterions 
there, and at the same time the repulsive intra-chain in- 
teractions are small. The overcharged branch is only ac- 
cessible if the chemical potential is further reduced at the 
same time that the screening increases. For most cases 
with sensible sizes for the counterions the reversal in di- 
rection of the parametrized curve comes always before 
the redissolution, and thus the redissolution is into the 
naturally charged state. 





I 

V Charged 












s^^* b 






r Neutral 








^^-''"''overcharged 
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FIG. 3. Diagram of fj, — n 2 showing the regions of expanded 
(gray) and collapsed (white) conformations. The solid curves 
give the relation of fi(n) using the Full Debye-Huckel model 
for different size of the ions a. The first line (a) is for a = 0, 
the second line (b) is for a = 0.2 and the third curve (c) is 
for pure entropy (i.e., neglecting the electrostatic interactions 
of the multivalent ions in the solution, and mathematically 
equivalent to a = oo). The three lines are indistinguishable 
at low concentrations (small values of k 2 ). 

It is interesting to note that at high enough salt con- 
centrations and high valences, the salt might re-associate 
on its own. In this case the trajectory in the /i — K 
plane has a different shape than that provided by the 
Debye-Huckel equation. Indeed, while the results from 
the Debye-Huckel model are know to provide good ap- 
proximations for the free energy of the salt in solution 
they do not properly describe the individual dissociated 
species. The upturn of the chemical potential essentially 
indicates a phase transition for the salt/solvent system. 
In reality there is no phase transition but there is a finite 



a fraction of ion pairs formed in the solvent (even if they 
do not reassociate chemically). We have not presented 
here the curve in the /i — k 2 diagram when the associa- 
tion of multivalent ions in the solution is present, as-is 
modeled, for example, in the classical Bjerrum modcllia, 
or more recently by Fisher, et al.EJ This association gives 
large corrections to the Debye-Huckel law (for the free 
energy of the individual species of ions) and in our case 
provides a /i — k curve that monotonically approaches a 
horizontal line (a limiting /i value), with increasing in- 
verse screening length. These effects are particularly im- 
portant for large valence salts as in our case and we-.will 
discuss them in more detail in a future publicatior£3. 

V. COMPARISON WITH EXPERIMENTAL 
INFORMATION 

With the previously obtained solutions we now con- 
trast our results with the currently known experimental 
situation. We consider the shape of the phase diagram, 
the precise location of the collapse transition and the pre- 
dicted structure of the polymers in both of the phases. 

It should be clear that we have already recovered the 
shape of the B-C transition. Indeed, the (approximate) 
equation that define the transition does not have any de- 
pendence on the monomer concentration, and therefore 
produces a flat line in the cf> — m, or <f> — jj, diagrams. 
The transition is defined by the choice of single chains to 
be in one of two states that must be in equilibrium with 
the surrounding environment. This environment, repre- 
sented here by means of an effective chemical potential 
and screening length, depends strongly on the amount 
of multivalent salt and only very slightly on the polymer 
amount. The more precise expressions for the chemi- 
cal potential recovers a small dependence in the polymer 
concentration. 

There are two important ways in which important 
changes on the qualitative form of the diagram can 
arise. The addition of large amounts of monovalent salt 
introduces new features in the phase diagram, not show 
in the scheme of Fig. I. In the introduction we only 
mentioned the effects on the A — B transition, but fur- 
ther addition of monovalent salt, creates a different envi- 
ronment of the chains, and changes the transition points 
and the shape of the transition curve. Secondly, it is clear 
that as we continue increasing the amount of polymer in 
the solution, the interactions between the chains start to 
become important and the transition (if any) is a much 
more complicated phenomena. From the point of view of 
our theory it is clear that the increase in concentration 
of monovalent salt and monomers bring about a break- 
down of the assumption that the chemical potential for 
the overcharging is dependent only on the amount of mul- 
tivalent salt. The effective environment of the chain be- 
comes more complicated and the chemical potential will 
now contain important terms coming form the concen- 
trations of all species and from many-body interactions. 



8 



Let us consider two concrete numerical examples. Ex- 
periments performed on polystyrene sulfonated at room 
temperature with N — 4 * 10 3 , Bjerrum number of 
B = 2.87 (b = 0.5 ran) and an effective ion radius 
a/b = 0.2, show a transition line at a concentration of 
to = 0.2M of LaCl^. With these parameters as input, 
our equation predicts redissolution at a salt concentra- 
tion of to = Q..1M. Secondly, experiments on double 
stranded DNAa with very large number of base pairs (we 
take N = 10 4 ) and B — 4.2 (b = .17 nm)in spermine 
(z = 4), a = 1. obtained redissolution at m = AM, 
while our equations predict a transition at m = 0.04M. 
This compares favorably with the experiments, especially 
if we take into account the rough approximations done 
in the evaluation of the energies of the system. 

According to the theory presented, in phase B, the 
polymer chains are almost neutral with a collapsed con- 
formation. In phase C, the chains are expanded and 
charged. This suggests, besides other techniques, to con- 
firm these predictions by means of scattering experiments 
that test the structure of the chains, and of osmotic pres- 
sure measurements that can determine the amount of free 
counterions in eajcJa-state (as it has been done for semi- 
dilute solutionstSl 3 !!) . 



VI. CONCLUSIONS 

The redissolution transition observed in multivalent in- 
duced precipitated polyelectrolytes with further addition 
of multivalent particles or salt was predicted extending 
a previously developed two state modelEl to deal with 
large salt concentrations. The electrostatic energies of 
both, collapsed and expanded-coiled, states were com- 
puted considering the finite size of the ions and monomers 
condensed along the chains, and using a mean field ap- 
proach for the non-condensed ions. We neglected the 
non-ionic short range interactions of the ions with the 
solvent, and assume zero monovalent salt. 

We found that the redissolution is determined by the 
properties of the ionic solution. We calculated the effec- 
tive charge of the chains for the different thermodynamic 
states of the chain. 
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